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We use the recently introduced U(N) framework for loop quantum gravity to study the dynamics 
of spin network states on the simplest class of graphs: two vertices linked with an arbitrary number 
N of edges. Such graphs represent two regions, in and out, separated by a boundary surface. We 
r*j . study the algebraic structure of the Hilbert space of spin networks from the U(N) perspective. In 

r-{ • particular, we describe the algebra of operators acting on that space and discuss their relation to the 

' standard holonomy operator of loop quantum gravity. Furthermore, we show that it is possible to 

, make the restriction to the isotropic/homogeneous sector of the model by imposing the invariance 

^s^j ■ under a global U(N) symmetry. We then propose a U(N) invariant Hamiltonian operator and study 

the induced dynamics. Finally, we explore the analogies between this model and loop quantum 
i— i, cosmology and sketch some possible generalizations of it. 
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Introduction 

Loop quantum gravity (LQG) presents a rigorous framework towards the quantization of general relativity. The 
main issues faced by this quantization scheme are the precise dynamics of the theory and the derivation of its semi- 
classical limit. There has been tremendous progress on these topics since the original formulation of the theory, but 
a punch line is still missing. Our main motivation for the present work is the study of the LQG dynamics. Since 
Thiemann's original proposal of a well-defined Hamiltonian constraint operator for LQG [l|, 0] , there has been various 
more recent proposals among which we point out the algebraic quantum gravity framework Q and the spinfoam 
approach e.g. 0, @- Our main source of inspiration for our current approach is the recent model introduced by 
Rovelli and Vidotto @. The logic behind their model is to implement the LQG dynamics on the simplest non-trivial 
class of spin network states, thus constructing a first order truncation of the full theory. They considered spin network 
states based on a fixed graph with two vertices related by four edges, so that their model can be called "tetrahedron 
LQG model" . Very interestingly, it was shown that this model can be understood as reproducing a cosmological setting 
in LQG and leads to a physical framework very similar to loop quantum cosmology [fl 0] ■ ^ was a l so shown that the 
same procedure can be successfully applied to the current spinfoam models @. This "dipole quantum cosmology" is 
the starting point of our work. 

We consider the generalization of the Rovelli- Vidotto model to spin network states based on a graph with still 2 
vertices but now with an arbitrary number N of edges. From their viewpoint, this should allow to introduce more 
anisotropy /inhomogencity in their model. Here, we start anew with a thorough study of the algebraic structure of the 
Hilbcrt space of spin networks on the 2- vertex graph. To this aim, we use the recently developed U(V) framework for 
SU(2)-intertwiners Considering our spin networks as two intertwiners (one for each vertex) glued by matching 

conditions representing the N edges between the two vertices, we use the operators of the U(V) formalism to act on 
our Hilbert space. From there, we introduce two improvements on this standard setting. First, we give the explicit 
link between the U(V) operators and the usual holonomy operators of LQG, which provides a dictionary to translate 
between the standard approach to LQG and the more recent U( N) framework. Second, we identify a U(V) symmetry 
generated by operators acting on the coupled system of the two vertices and which reduces the full space of arbitrary 
spin network states to the reduced space of homogeneous/isotropic states. 

Then we finally introduce U(V)-invariant operators, which respect the isotropy /homogeneity of the quantum states, 
and use them to build a U(V)-invariant Hamiltonian operator encoding the dynamics of our 2- vertex model. We study 
the algebra of our operators and propose an interpretation of our creation operator as a "black hole" creation operator. 
Finally, we study the spectral properties of our Hamiltonian and discuss in detail the interplay between our framework 
and loop quantum cosmology. We conclude with the limitations of our model and its possible generalizations in order 
to take into account more physical situations. 

I. THE U(N) FRAMEWORK FOR LQG INTERTWINERS 

Let us focus on the intertwiner spaces for the Lie group SU(2) and define the Hilbert space of intertwiners between 
N irreducible representations of spin j±, .., j'at : 

It j N = lnv[V jl ® .. g> V jN ]. (1) 

The object that will be most relevant to us is actually the space of intertwiners with N legs and fixed total area 
J = E* 3i ■ 

We also introduce the full Hilbert space of iV-valent intertwiners: 

V/v 0//,......., 0V/y. (3) 

The key starting point to our dynamical model is that there is a natural action of U(V) on the intertwiner space 

Hn [1], that more precisely the intertwiner spaces Ti^ carry irreducible representations of U(iV) [l(| and finally that 
the full space 7-Ln can be endowed with a Fock space structure with creation and annihilation operators compatible 
with the U(iV) action p]]. 
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This U(iV) formalism is based on the Schwingcr representation of the su{2) Lie algebra in terms of harmonic 
oscillators. Let us introduce 2N oscillators with creation operators <ii, hi with i running from 1 to N: 

[oi, a]] = [bi, &j] = Sij , [en, bj] = 0. 

The generators of the SU(2) transformations acting on each leg of the intertwiner are realized as quadratic operators 
in terms of the oscillators: 

Jf = I(oJai - b\h), Jt=a\bi, Jr=a t bl E i = (a\a i + b\b i ). (4) 

The Jj's satisfy the standard commutation algebra while the total energy Ei is a Casimir operator: 

[J?,jf] = ±J±, [jf,Jr]=2J*, [E i ,T i ]=0. (5) 

The correspondence with the standard \ j, m) basis of su(2) representations is simple. The eigenvalue m, of Jf is given 
by the half-difference of the energies between the two oscillators, while the total energy Ei gives twice the spin, 2ji, 
living on the i-th leg of the intertwiner: 

\n a ,nb)oH = \\i n a + nb), ^(n -ra&)), \j,m) = \j + m,j - m) H (6) 
Intertwiner states are by definition invariant under the global SU(2) action, generated by: 

j z = J2 j >> J± = J2 J t- ( 7 ) 

i— 1 i 

Then, operators acting on the intertwiner space need to commute with these operators too. The simplest family of 
invariant operators was identified in Q and are quadratic operators acting on couples of legs: 

En = a\ aj + b\bj, El = E ji . (8) 

The starting result is that these operators are invariant under global SU(2) transformations and form a u(N) algebra: 

[J, Ekj] = 0, [E tj ,E kl ] = 6 jk E u - S u E kj . (9) 

The diagonal operators B, = En form the Cartan sub-algebra of u(N) , while the off-diagonal operators Eij with i ^ j 
are the raising and lowering operators. As said earlier, the generators Ei give twice the spin 2ji while the U(l) Casimir 
E = J^i Ei w iU gi ye twice the same area, 2 J = ^ 2ji. Then all operators E^ commute with the U(l) Casimir, thus 
leaving the total area J invariant: 

[E ij ,E] = 0, (10) 
The usual SU(2) Casimir operators have a simple expression in terms of these u(N) generators 1 : 

{■h? = \Ei ( f + x ) - V ^ i. ^ ■ Jo) = \ E >jE ]t - ~E % E, - ±Ei. (11) 

Moreover, the explicit definition of the E^j's in terms of harmonic oscillators leads to quadratic constraints on these 
operators as shown in [Io| . Here, we give a slight generalization of these constraints: 

where we have assumed that the global SU(2) generators J vanish. This can be interpreted as the square of the matrix 
E^ giving back itself up to a factor which is its trace E divided by 2 up to an ordering factor (N — 2). Classically, 



1 Let us point out that the case % = j of (J; • Jj) does not give back exactly the formula for (J;) 2 due to the ordering of the oscillator 
operators. The two formula agree on the leading order quadratic in Ei but disagree on the correction linear in Ei. 
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considering a (Hermitian) matrix M satisfying the similar constraint M 2 = |(TrM)M, it is straightforward to show 
that M has only a single eigenvalue with degeneracy 2. We will see below that this interpretation actually holds at the 
quantum level. This viewpoint of seeing the u(N) framework for intertwiners as the quantization of a classical matrix 
model will be investigated in detail elsewhere [12| • These quadratic constraints on the £7y operators lead to non-trivial 
restrictions on the representations of u(N) obtained from our construction. To solve them, the method used in [Io| is 
to apply them to a highest weight vector. Thus let us introduce a highest weight vector v, such that it diagonalizes 
the generators of the Cartan sub-algebra E{V — Uv with the eigenvalues h > I2 > h > •• > and that it vanishes 
under the action of the raising operators Fy v = for all i < j. Applying the quadratic constraints given above to 
such vector v, one shows that the eigenvalues must necessarily be of the following form, [h,h, h~, In] = [h 1, 0, .., 0]. 

Coming back to the intertwiner spaces, we get that the intertwiner space Jiff at fixed total area J = J^i 3i carries 
an irreducible representation of U(iV) with highest weight [J, J, 0, .., 0] with I = J, as was shown in This highest 
weight vector corresponds to a bivalent intertwiner with jx = j% = ^ and 33 = ■■ = Jn = 0. Then one acts on it with 
u(JV) operators to get all possible intertwiner states with the same total area J. Finally, the value of the quadratic 
V(N) Casimir is easily obtained from the previous constraints: 

E E ik E u = E + N - 2) = 2 J( J + N - 2). (13) 

i,k ^ ' 

And we can also compute the dimension of in terms of binomial coefficients using the standard formula for U(iV) 
representations : 

n -a- 1 ( N+J-l\ { N+J-2\ (N + J-l)l(N + J-2)\ 

D N ,j = dim7&' = — ^ j j j j = J!(J+1)! ( jV _ 1)!(jV _ 2) , ■ (14) 

Next, we introduce annihilation and creation operators to move between the spaces T-L^ with different areas J 
within the bigger Hilbert space of all intertwiners with N legs [ll[. We define the new operators: 

Fij = (cubj - cijbi), Fji = -Fij. (15) 

These are still invariant under global SU(2) transformations, but they do not commute anymore with the total area 
operator E. They nevertheless form a closed algebra together with the operators E^: 

[Eij,Eki] = SjkEu — SuEkj 

[E^Fu] = S a F jk - SikFjt, [Ei^Fl] = 5 jk Fl - S Jt F^ k , (16) 

[Fij,F kl ] = SikEij — SuEkj — SjkEu + SjiE^i + 2(SikSji — SuSjk), 

[Fiji Em] = 0, [Ft, Ft] =0. 

The annihilation operators F^ allow to go from Ti^ to while the creation operators F^ raise the area and 

go from t° ^7v +1 '- ^ n —I' these operators were mainly used in order to construct coherent states transforming 



consistently under V(N) transformations. Here, we will focus on their use as building blocks for the dynamics of 
intertwiners. 

These new operators satisfy new quadratic constraints together with the E^ . Starting with their explicit definitions 
in terms of oscillators, it is easy to get on the intertwiner space (assuming global SU(2) invariance, J = 0): 

E4^=4f> E^4 = 4(f+^- 1 )- (173 
E E * F * = F * (| - x ) - E F * E v = F « (f + N 2 ) ' ( 18 ) 

E F ikFjk - Eij ( f + l) , E F * F l = + 25 «) ( l + N - 1 V ( 19 ) 



k 



These relations will be useful later for the analysis of the dynamics of our 2- vertex toy model. It is also interesting 
to look at them from a classical matrix perspective. Indeed ignoring the "quantum" ordering terms, we consider 
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the classical Hcrmitian matrix AI = M' corresponding to the operators E^ and the classical antisymmetric matrix 
Q = — Q* corresponding to the operators and the quadratic constraints become simply: 

A/ 2 = AM, QQ = -XM, 

QM = XQ, MQ = XQ, (20) 

where the eigenvalue is defined in terms of M as A = (TrM)/2. These conditions correspond exactly to the relations 
between the matrices p and z introduced in in order to define U(iV) coherent states. This strongly suggests that 
this U(iV) framework can be derived as the quantization of a matrix model. This will be investigated in detail in 
future work [l2j]. 

II. THE 2- VERTEX GRAPH AND THE U(JV) SYMMETRY 
A. The 2- Vertex Algebraic Structure 

Let us consider the simplest non-trivial graph for spin network states in Loop Quantum Gravity: a graph with two 
vertices linked by N edges, as shown in ngffl This is a generalization of the simplest tetrahedral model introduced 
by Rovelli and Vidotto in Q and which was shown to be related to models of quantum cosmology Q • Their model is 
made of two vertices related by exactly four edges, which can be interpreted as two tctrahedra glued together. Here 
we consider the more general case of an arbitrary number N of edges. 




FIG. 1: The 2-vertex graph: the two vertices a and /3 are linked by N edges. 

Naming the two vertices a and /3, we now have two intertwiner spaces attached respectively to a and (3 with their 
own operators e[^\f^ and E^' ,F^P . The total Hilbert space of these uncoupled intertwiners is the tensor product 
of the two spaces of A'-valent intertwiners: 

U m = U N ®H N = H^^H^ = U j? ,.., j% 8%,..^. (21) 

This space carries two decoupled U(iV)-actions, one acting on the intertwiner space associated to the vertex a and 
the other acting on f3. 

However, considering the 2-vertex graph, the two intertwiner spaces are not independent. There are matching 
conditions [ljj imposing that each edge carries a unique SU(2) representation, thus the spin on that edge must be the 
same seen from a or (3 i.e. jf = jf . This translates into the fact that the oscillator energy for a on the leg i must be 
equal to the energy for j3 on its i-th leg: 

Si = E\ a) - E^ = 0. (22) 

This set of conditions is stronger than requiring that the total area of a is the same as /?, even though this is a 
necessary condition. Then, the Hilbert space of spin network states on the 2-vertex graph is much smaller than the 
decoupled Hilbert space W®2- 

2 n 0/c:.. M (23) 

Then, in order to define consistent operators acting on 2 "H, we not only have to check that they are invariant under 
global SU(2) transformations, but we also have to check that they commute (at least weakly) with the matching 
conditions Si. 
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We can construct operators deforming consistently the boundary between a and j3. We introduce the following 
symmetric operators: 

We check that they commute with the matching conditions 

Vi,j,k, [£k,eij] = [£k,fij] = 0. (25) 

Calling £ = J2 t E[ a) = J2i s f ) the operator giving (twice) the total boundary area on our Hilbcrt space 2 H satisfying 
the matching conditions, the operators e^'s preserve the boundary area while the 's will modify it: 

[E, eij }=0, [B,/i i ] = -2/i J) [£?,/t ] = +2/ t. (26) 

More precisely, the operator increases the spin on the i-th edge by +^ and decreases the spin of the j-th edge. 
The operator /„ decreases both spins on the edges i and j, while its adjoint fjj increases both spins by \. These 
operators generate the deformations of the boundary surface, consistently with both the SU(2) gauge invariance and 
the matching conditions imposed by the graph combinatorial structure. They will be the building blocks for the 
dynamics of spin network states on the 2-vertex graph. 

Nevertheless these operators e.%j,fij,flj do not have simple commutators with each other. Moreover, they are not 
enough to generate the whole Hilbert space 2/ H of spin network states from the vacuum state |0). Indeed, they arc 
symmetric in a O f3 and will only generate states symmetric under the exchange of the two vertices. To explore 
the whole Hilbert space 2 H and generate the whole space of gauge invariant operators on the 2-vertex graph, we 
also need operators that act on a single vertex and change states without affecting neither the representations on the 
edges, nor the intertwiner state living at the other vertex. Natural candidates for such operators acting on j N 

are the Eg>E£\ which are directly related to the usual scalar product operators ■ by (jTTj) . They change 
the intertwiner living at the vertex a without changing the intertwiner living at j3. Actually, it is straightforward 
to see that combining such local operators E\ a \ ' E^f' and E^\E^'E^' with the coupled symmetric operators 
e ij 5 fij i fjj allow to go between any two states in the Hilbert space 2 H of spin network states on our 2-vertcx graph 
and thus generate all gauge invariant operators (i.e the algebra of holonomy and grasping operators) on that graph. 



B. The U(JV) Symmetry, the Invariant Subspace and Homogeneous States 

As we have said above, the matching conditions £k break the U(7V)-actions on both vertices a and (3. Nevertheless, 
we can see that the £ k generate a U(l) symmetry and that they are part of a larger U(7V) symmetry algebra. Indeed, 
we introduce the operators: 

£ l3 = E^-Ef =E\f> -(i$>)t. (27) 
It is straightforward to compute their commutation relations and check that these operators form a u(iV) algebra: 

[£ij,£ki] = 5jk£a — 5u£kj- (28) 

The diagonal operators are exactly the matching conditions £kk = £k an d generate the Cartan abelian sub-algebra of 
u(N). 

These £V,'s generate U(A) transformations 2 on the two intertwiner system that act in as (U,U) with the 

transformation on ft being the complex conjugate of the transformation on a. In the following, unless otherwise said, 
we will refer to the transformations generated by the £y simply as the U(iV) action for our 2-vertex setting. 

Looking closer to the interplay between the matching condition and these new u(N) generators, 

[£ij,£k] = Sjk£ik — 5ik£kj = {&jk ~ Sik)£ij, (29) 

two points are obvious: 



2 There is another set of U(iV) transformations acting on the coupled system as (U, U) and generated by: 

E M = E [a) + E W . 
ij ij ij 

They satisfy the same commutation relations, [EY^',Ej^^] = Sj^E^' — S^E^' , but they will not be relevant for our construction. 
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• The operators £{j are not fully compatible with the matching conditions and they do not act on the 2-vertex 
Hilbcrt space 2 H. Thus they do not generate a non-trivial U(AT)-action on 2 H. Indeed, applying the operator 
Sij on a vector satisfying the matching conditions, £k IV') = 0, we obtain a vector which does not satisfy the 
matching conditions unless it is trivial: 

£fc£sj |</>) = [£k,£ij] \i>) = (Sik - $jk)£ij \ip}- 

• We can nevertheless look for vectors in 2 T-L which are invariant under this U(iV) action, \ijj) = for all 
In particular, they will satisfy the matching conditions (given by the case i = j). 

Following this line of thought, we introduce the subspace of spin network states which are invariant under this 
U(iV)-action: 

The first equality holds because the matching conditions £k are part of the u(iV)-algebra, while the second equality 

is simply the definition of the full decoupled Hilbert space H®2- Now taking into account that the spaces T-L^ are 
irreducible U(iV)-represcntations, requiring U(-/V)-invariancc imposes that the two representations for the two vertices 

are the same, J a = Jp, but furthermore there exists a unique invariant vector in the tensor product Hffi <8> n}^ . We 
will call this unique invariant vector | J) and we will construct it explicitly in terms of ey and fij in the section 1111 Bl 
below. 

What's important here is that, by imposing U(7V)-invariance on our 2-vertex system, we obtain a single state | J) 
for every total boundary area J: 

2 H mv = ($C\J). (31) 

Jen 

This means that we have restricted our system to isotropic states, which are not sensitive to area-preserving defor- 
mations of the boundary between a and (3. Let us point out that these U(7V)-invariant states are both isotropic and 
homogeneous. Indeed, they are isotropic since all directions are equivalent and that the state only depends on the 
total boundary area, and they are homogeneous because the quantum state is the same at both vertices a and j3 of 
the graph. This has two interesting consequences: 

• It allows to realize the reduction to the isotropic/ homogeneous subspace used in the tctrahcdral model by Rovclli 
and Vidotto 0,0] to separate the homogeneous dynamics of the global boundary area J from the inhomogeneous 
degrees of freedom corresponding to the area-preserving variations of the individual spins j, living on its edge. 
Here, we've shown that the projection on isotropic/homogeneous states is achieved by a straightforward XJ(N)- 
group averaging. 

• This logic can be applied more generally to loop quantum cosmology, which is based on a symmetry reduction at 
the classical level and a quantization a la loop of this reduced phase space. A challenge is to derive loop quantum 
cosmology from the more general realm of loop quantum gravity. A first important step in this direction was 
achieved by the tetrahedral model of Rovelli and collaborators where they could derive a simple isotropic 
and homogeneous quantum cosmology from loop quantum gravity on the double-tetrahedron graph i.e the 2- 
vertex graph with N = 4 edges. Here, we would like to point out that the present U(7V) symmetry allows 
to perform the symmetry reduction to homogeneous (isotropic) states directly at the quantum level. It seems 
that this has been an important missing ingredient in order to understand the symmetry reduction from loop 
quantum gravity to loop quantum cosmology. 

In the following sections, we will impose this U(iV)-invariance, work on the invariant subspace of homogeneous 
states and look for U(iV)-invariant dynamics which would act on this invariant space. In particular, such dynamics 
will automatically be consistent with the matching conditions. 

C. Back to Holonomy Operators 

In order to understand the precise relation between the new U(iV) formalism that we use here and the more standard 
framework of loop quantum gravity, it is interesting to investigate the link between our operators and fij and the 
usual holonomy and derivation operators of loop quantum gravity. 



n 



N 



J a ,J B 



(30) 
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To define the action of the holonomy operator, let us quickly review the definition of the spin network functionals. 
We consider the 2- vertex graph and choose to orient all the edges in the same way from the vertex a to the vertex /?. 
The wave-functions are functions of N group elements in SU(2) which satisfy a SU(2)-invariance at both vertices: 

^(9i,-,9n) =ip(h~ 1 g 1 hp,..,h~ 1 g N hp), Vh ai hp g SU(2) x2 . 

Then the (simplest) holonomy operator acts on a couple of edges (ij) by multiplication: 

{x m 1>)(gi,~,9N) = Xi{9i9j l )^{9r,..,9N), (32) 

where X| is the spin-i character, which takes the trace of the SU(2) group element in the fundamental representation 

of spin This operator is clearly Hermitian and symmetric under the exchange i j. 

To understand the action of the holonomy operator, we need the Clcbsh-Gordan coefficients corresponding to 
tensoring an arbitrary spin j with the spin |. Let us then consider the tensor product V 3 : <g> V? which decomposes 

as is well-known as y^'+s © V^~^. Let us write J and j for the su(2) generators acting respectively on V 3 and Vi 
while we call J = J + j the generators of the su(2) transformations acting on the coupled system. The scalar product 
operator J- j has two eigenvalues, | and — t » which correspond to the two proper subspaces and V^~^ . We 

use the oscillator basis where we drop the index OH, 

\n,2j-n)® |e,l-e), 

where n running from to 2j is the energy of the a-oscillator in V-* while e = 0, 1 labels the basis of Vi. We remind 
that the spin j is given by the half-sum of the energies of the two oscillators and the magnetic momentum m is given 
by half- difference of these energies, so that this state corresponds to \j, m = n — j) ® ||, e — |) in the standard SU(2) 
basis. We can compute explicitly the action of the scalar product operator in this oscillator basis: 

J-J\n,2j-n)®\l,0) = i(n - j)\n, 2j - n) <g> |1, 0) + 1 y/(n + l)(2j - n) \n + 1, 2j - n - 1) <g> |0, 1), (33) 
J-)\n,2j-n) ® |0,1) = -(j - n)|n, 2j - n) <g> |0, 1) + -y 7 ^ - n + 1) \n - 1, 2j - n + 1) ® |1, 0). (34) 

From there, it is easy to check that the subspace V^+l is spanned by the following basis vectors which are eigenvectors 
of J ■ ) for the eigenvalue + 1 : 

\n,2j + l-n) = (v 7 ^ |n - 1, 2j - n + 1) ® |1, 0) + ^/2j + l-n \n, 2j - n) ® |0, 1)) , (35) 

for which the eigenvalue of J z is m = (n — j) — \ . Here the energy level n runs from to 2j + 1. On the other hand, 
the basis of the subspace V^~^ is given by the states orthogonal to the previous ones: 

\n-l,2j-n) = (y/2j +T^|n - 1, 2j - n + 1) ® |1, 0) - yfc \n, 2j - n) ® |0, 1)) , (36) 

for which the eigenvalue of J z is still m = (n — j) — h but now the energy level n runs from 1 to 2j. We also give the 
reverse mapping: 

|n,2j-n><8>|0,l) = ^= (V^fTl^ n \n, 2j + 1 - n) - yfti \n - 1, 2j - n>) , (37) 

|n,2j-n)®|l,0) = (Vn + 1 \n + 1, 2j - n) + ^/2j - 7i \n, 2j - 1 - n)) ■ (38) 

To express the holonomy operator in terms of the u(iV) operators, it is more convenient to express the previous 
decomposition formula in terms of the oscillator operators: 

\n, 2j - n) ® |0, 1) = ^JTT ^ ~ ^ |n ' 2j ~ n) ' (39) 
\n, 2j - n) ® |1, 0) = (a 1 + &) I", 2j - n). (40) 
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FIG. 2: Holonomy operator x lJ acting on the couple of edges (ij), which sends the spins to the four possibilities 

(ji +Vj) with rji,r]j — To study the precise action of this holonomy operator, we focus on its action Xa an d Xp on 

the individual vertices a and /3. Then we glue back the two resulting intertwiners while respecting the matching conditions. 



Now coming back to the holonomy operator, it can be seen to act independently on each intertwiner at the vertices 
a and (3, as shown on figiH Focusing on the state at the vertex a, the holonomy operator acts on the legs i and j of 
the intertwiners by tensoring the state on these legs with the singlet state (|1, 0) (g> |0, 1) — |0, 1) ® |1, 0)) of the spin-i 
representation. This corresponds in figure [2] to the insertion of the intermediate edge labeled by the spin-i between 
the legs i and j. Dropping the index a, this action reads : 



X ( a j) ■ ®«k — ► (vi®Vj® (|1,0)® |0, 1> — |0, 1> <g> |1,0»^ (g) v k 

fe=l ^ ' 

= ((vi ® [1, 0>) ® (vj ® |0, 1)) - (yi ® |0, 1)) ® (yj <g> |l,0))j 0u fe . (41) 

The crucial point is that this map is SU(2)-invariant since the singlet state is invariant under SU(2). Thus it is a map 
between intertwiner spaces: 

X W> : Inv[V Jt ® Vfc (g) ^] — ► Inv[(^ ® Vj) ® ® Vj) (g) V},]. 

ft fc 

This map decomposes into four components following the decomposition into irreducible representations of the tensor 
products on the edges i and j : 

(V Jt ® Vj) ® ( V}, ® Vj ) = (V* ® v., +% ) . 

Vi , Vj =±l/2 

This decomposition is done using the previous Clebsh-Gordan formulas (|39H40[) . The holonomy operator acts like this 
on both intertwiners on the vertices a and /3 and then glues them by requiring that the edges i and j carry the same 
spins ji + rji and jj + rjj seen from both vertices a and /3 following the matching conditions. Putting all the pieces 
together, we obtain a very simple expression for the holonomy operator in terms of the oscillator operators acting on 
both a and /3, which can be written in terms of eij and operators: 



1 

The four terms correspond to the four sectors of the holonomy operators, while the inverse-squareroot factors come 
from the l/ v / 2j + 1 factors of the Clebsh-Gordan coefficients. Expressed as such, the holonomy operator is clearly 
SU(2) invariant, Hermitian and consistent with the matching conditions. Finally, the operators (Ek + 1) are Hermitian 
and strictly positive, thus we can legitimately consider their squareroot and inverse. 

This formula is surprisingly simple and expresses the holonomy operator from loop quantum gravity in terms of the 
basic operators of the U(iV) formalism. It establishes a dictionnary between the usual formulation of loop quantum 
gravity in terms of holonomy and derivation (grasping) operators and our formulation in terms of the u(iV) operators 
and the creation/annihilation operators Fij,F^. 

Next, we check the previous formula on bivalent intertwiners, which are the simplest states of the 2-vertex Hilbcrt 
space of spin networks. A bivalent intertwiner state is given by setting the spins to ^ on two chosen edges i and j and 
to on all other edges. This defines a unique spin network state given by the tensor product of the same singlet state 
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of the spin--^ representation put on both vertices a and (3. It is also obtained by applying the holonomy operator in 
the spin--^ representation on the loop (ij) formed by the couple of edges. Ignoring all the other edges k ^ i,j, the 
explicit expression for this state reads: 



K4 



.7 



^7 ( E^ 1 )"!"' J — n)®\J — n, n)\ ® fe(-l)*\p 

\n=0 / \p=0 



J-p)® \ J-p,p) 



(43) 



where the first vector in both a and j3 terms gives the state living on the edge i while the second vector defines the 
state on the edge j. This state is normalized, (xj.\x±) = 1- Moreover, this is a straightforward tensor product state. 

Thus to derive the action of the operators e$j = E^'Ey and /y = F^F^' on this state, we can drop the index a 



(and 0) and look at the action of the operators Ey and Fij on the singlet state ^n=o( — l)™! 71 ; J — n) ® \ J — n,n) . 
Going back to their expressions in terms of oscillator operators, we easily compute: 



Eij V](-l)"|n, J -n)®\J — n,n) 

n=0 
J 

Eji ^(-l) n |n, J -n)®\J -n,n) 

J 

Fij ^(-l) n |n, J-n) ® \J-n,n) 

./ 

4 £(-!)>, </-n)® \J~n,n) 



(44) 



j-i 



-(J + 1) 5^(-l) n |n, (J - 1) - n) ® |(J- 1) - n, n) 

J+l 



n=Q 



Therefore, in the sum (/^ + + e 3 i + fij) giving the holonomy operator, only the two terms fij and changing 
the total area boundary contribute while the two other operator and eji vanish on the bivalent intertwiner state. 
This gives: 



X 



(ij) 



1 



1 



VE~TTy/E~+T 



X4 



Xi^l ) + 



Xi+i 



(45) 



And we recover the standard action of the holonomy operator on a single loop state. Let us point out that this action 
allows to easily derive the eigenvalues and eigenvectors of the holonomy operator. Indeed, the eigenvectors are given 
by states whose coefficients are themselves the evaluation of the characters on group elements with arbitrary but fixed 
class angle 9 : 



J6N 



JEN 



(46) 



where the evaluation of the characters are: 



sin(J+ 1) 



X o(0) = 1, X i{0) = 2cos0. 



As we can see, the spectrum of the (one-loop) holonomy operator is bounded in absolute value by 2. Moreover, the 
physical interpretation of these eigenvectors is straightforward. In the limit case 8 — > 0, the spin network functional 
corresponding to the state ^2j(J + 1) \Xj/2) is simply the distribution S(gig~ 1 ) defining the flat connection state in 
loop quantum gravity. More generally, for an arbitrary angle 6, the corresponding spin network functional is the 
distribution fixing the class angle of the holonomy gig J 1 to 8. 

To conclude this section, we would like to give the reverse formula of the equation (|4"2"|) giving the expression of 
the holonomy operator in terms of the u(N) operators. To recover the individual components fij, f}j from the 
holonomy operator x ■, we need to project on each of the four sectors of the holonomy operators. This can be done by 
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combining with specific insertions of Ei and Ej. For instance, we can select the transition (ji,jj) — > + ij + i) 
by the following combination : 



i ([*, ] + [*, •] + ft, ■] + 1) *<«> = 4 



1 ; . (47) 

y/Ei + ly/Ej + 1 

We can adapt this formula to select the other transitions and select the other components of the holonomy operator: 

+ = . 



This concludes the dictionary between the holonomy operators and the operators of the V(N) formalism. 

Now that we have explored in detail the algebraic structure of the space of spin network states on the 2-vertex 
graph, we will look at the U(iV)-invariant operators acting in our Hilbert space 2 T-L and will define U(-/V)-invariant 
dynamics that evolve homogeneous states. This will allow us to make a link between our formalism and the standard 
framework of loop quantum cosmology. 



III. DYNAMICS ON THE 2-VERTEX GRAPH 
A. The algebra of U(iV)-invariant operators 

In the previous section, we have introduced the U(iV)-invariant space 2 T-Li nv of homogeneous states spanned by basis 
vectors | J) labeled by solely the total boundary area. Now, we would like to study the structure of U(7V)-invariant 
operators that will act on this space 2 Hi nv of invariant states. 

The most obvious U(7V)-invariant operator is the total boundary area operator E itself. It is defined as E = E^ a ' = 
on the space % of spin network states satisfying the matching condition. It is direct to check that it commutes 
with the u(N) generators: 

[£ij,EV>] = \e\?,eW] = = -\Ef,E^\ = [E ih E^]. (49) 
This total area operator is clearly diagonal in the basis | J), 

E\J) = 2J\J). 

Now, we need operators that could create dynamics on the space 2 Hi nv by inducing transitions between states with 
different areas. To this purpose, we use the operators £ij-,fij and introduce the unique linear combinations that are 
U(iV)-invariant: 

ij ij ij ij 

They obviously commute with the matching conditions since each operator and fij does. Moreover, it is straight- 
forward to show that they are also invariant under U(iV) transformations: 

[£ij,e} = 0, [£ij,f] = 0, [£ij,p} = 0. (51) 

These three operators e, /, p seem to be the only U(-/V)-invariant operators which are quadratic in E and/or F and 
satisfy the matching conditions (the action of the operators on both vertices has to be the same in order to preserve 
these conditions, so it is necessary for such quadratic operators not to mix the operators Eij, Fij and F^ in the same 
term). 

Using the quadratic constraints (|17H19j) satisfied by the operators Eij and Fij, we can show that e, / and p form 
a simple algebra: 

[e,f] = -2f(E + N-3) = -2(E + N-l)f, (52) 
[e,/t] = 2p(E + N-l), (53) 
[/,/t] = A(E + N)(e + 2(E + N -1)). (54) 
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Looking at these commutators, it seems natural to introduce a shifted operator e = e + 2(E + N — 1). Then the 
algebra reads: 

[ej] = -2(E + N + l)f, (55) 
[e,/t] = 2p(E + N + l), 
[/,/t] = 4(E + N)e. 

Written as such, it resembles to a SI2 Lie algebra up to the factors in E, which is an operator and not a constant. 



This is very similar to the 5I2 algebra Ep a \ F^ a \ F^ 1 ^ defined in [ll| and used to build the U(iV) coherent states 3 
Since /, f* play the role of lowering and raising operators in this algebra, we can use f* as a creation operator. Thus 
we introduce the states: 

\J)un^P J \o) = (t /•;;' /;; j |o>, (56) 

where the index un stands for un- normalized (or for U^) following the reader's preference). Since both the operator 
P and the vacuum state |0) are U(A r )-invariant, it is clear that the states \ J) un are also invariant under the XJ(N)- 
action. Moreover, it is easy to check that they are eigenvectors of the total area operator: 

F |</)un — 2(/ |«/)ttn; 

so that they provide a basis for our Hilbert space 2 T-Li n v of homogeneous states. Then, using the commutators (|52ti54|) . 
we can compute the action of the operators e, / and f> on these states 4 : 

e|J)„„ = 2J(N + J-2)\J) un , (57) 

e\J) un = 2{J+l)(N + J -l)\J) un , 

/I Am = AJ{J+l){N + J-l){N + J-2)\J-l) un , 

P\J)un = |J+l)u„. 

As expected, the states \ J) un diagonalizc e, while /' and / act respectively as creation and annihilation operators. 

The action of e provides a consistency check. Indeed, taking into account the u(N) constraints E^" — E^jp = 
imposed on our invariant space n.i nv , the operator e is simply equal to the quadratic U(iV)-Casimir on that Hilbert 
space : 

e\J) un = J2 E £ )E iP i J )«« = E4 Q H Q) l J )u„ = E(a) (^r+ N ~ 2 ) l J )™ = 2J(J + N-2)\J) un , (58) 



=,3 



In the previous work [ill , operators E p a \ , F^^ were constructed by contracting the operators E^\f>"\f^"^ with the 
Hermitian matrix pij and the antisymmetric matrix both defined in terms of the set of N given spinors {(z® , z})}j=i..jv : 



Up =2^ Pn ' Pij = z i z j + z i z j > 

ij 

F ( a )i = ^ z ..j?Wt z ..= 2 o 2 i_ 2 oi 

z — l 3 ij > l J — z i z j z ] z % ■ 

ij 

Then the three operators e\"\ i*^ ■ , form a sl2 algebra. In our framework, working on the two- vertex graph, we have replaced 

these two classical matrices p and z by the operators acting on the /3- vertex, pij —> and z;j — > F^y^. This basically amounts to 

replace the spinor labels by the annihilation operators of the /3-vertex, (z®,zj) — ► From this point of view, our 2-vertex 

framework generalizes the algebraic structure of the single intertwiner case and it is natural to find similar algebras. 
4 We use the commutation relations of e and / with /t to compute their action: 

J-l J-l 
e| J)un = efi J \0) m = ^/ tJ - 1 - fc [e,/t]/t'=|0) = 2 £(2fc + N - 1) /t J \0) = 2J(N + J - 2) /t », 



fe=0 



J-l 



/|J> U n = // tJ |0> = 53/ tJ ~ 1-fc [/./ t ]/ tk l°> = 4^(2fc + 7V)(2fe(fc + 7V-2) + 2(2fc + 7V-l)))/tJ-l|0>. 
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where the value of the Casimir operator was given earlier in (|13[) . Similarly, we can express the shifted operator e as 
a polynomial in the total area operator E on the invariant space 2 Tii n v 



e |V> = \{E + 2){E + 2(N - 1)) |t/>), W> e 2 U mv . (59) 

Using this action, we compute the norm of the states \J) U n by recursion: 

un(J\J)un = (Ol/ 7 - 1 f\J) = 4J( J + l)(N + J - 1)(N + J - 2) un (J - 1| J - 1)„„, (60) 
which leads us to the scalar product: 



(61) 



in terms of the dimension of the intertwiner space nffi given in (|14[) . Thus we can define normalized basis states: 

|J)s i^(7Tk^7 IJ> ""' (62) 

on which the action of all e, /, y operators is always quadratic in J: 

e|J) = 2(J + l)(iV + J- 1)| J), (63) 
/| J) = 2^J(J + 1)(N + J - l)(N + J - 2) | J - 1), 
/t|J) = 2 V( J + W + 2)(iV + J)(N + J - 1) | J + 1). 

From here, we see that it is possible to introduce renormalized operators that truly form a SI2 algebra. We define 
the new operators: 

Z = 



y/E + 2(N-l) ' V /E + 2(N-1) ' 

^- = 1 — f 1 — , (64) 

y/E + 2{N - 1) ^E + 2(iV-l) 

x = 1 f t 1 

+ ^ + 2(^-1) V / E + 2{N- 1) 

First, the inverse square-root is well-defined since i? + 2(7V — 1) is Hcrmitian and strictly positive as soon as N > 2. 
Moreover, these operators are still U(./V)-invariant since E is invariant too. We also have the Hermiticity relations, 
Z* = Z and XI=X+. It is fairly easy to compute the action of these renormalized operators on our | J) basis states: 

Z\J) = (J+1)\J), (65) 
X_|J) = y/j(J+l)\J-l), 
X+\J) = v/(J + l)(J + 2)|J+l). 

These operators clearly form a 5/(2, R) Lie algebra: 

[Z, X±] = ±X±, X_] = -2Z. (66) 

Furthermore, the action (|65p corresponds to the irreducible unitary representation of sl(2,M) with vanishing Casimir 
Q = — Z 2 + {X-X+ + X+X-)/2 = 0. This irrep belongs to the positive discrete series of irreps of sZ(2,R), which 
have a spectrum bounded from below for the operator Z (the interested reader can find more details on the unitary 
representations of sZ(2,R) in [l3jV 
This shows two important points: 

• The algebraic structure of the U(-/V)-invariant space 2 Hinv of homogeneous states is clear: it forms an irre- 
ducible unitary representation of sZ(2,R). The basis vectors | J) can be obtained by iterating the action of the 
creation/raising operator /' (or X + ) on the vacuum state |0). 
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• This algebraic structure does not depend at all on the number of edges N. Therefore, while working on 
homogeneous states, N might have a physical meaning but it is not a relevant parameter mathematically. On 
the other hand, we expect it to become highly relevant when leaving the U(iV)-invariant subspace and studying 
inhomogeneities. 

Since the algebra of the operators Z, X± is much simpler than the algebra of the operators e, f, f' , it is also convenient 
to give the reverse formula and write e, /, f' in terms of the s?2 operators. Using that Z is simply related to the total 
area E on the invariant space 2 T-Li n v, 

Z= — + l, E + 2(N-l) = 2(Z + N-2), 

we have the inverse relations, which hold on the Hilbcrt space 2 'Hi nv : 

e = 2Z(Z + N - 2) , (67) 
/ = 2 VZ + N -2X-^Z + N - 2, 

f = 2v*Z + N-2X+VZ + N-2. (68) 



Finally, to conclude this section, we introduce some further renormalized operators: 

i?4=i, -4/4, 4/ f i (69) 

Ve ve ve ve ve ve 

These operators obviously still commute with the U(iV) action. They are also well-defined since e is Hermitian and 
strictly positive. Their action on the basis states is the simplest, since they generate straightforward translations 
between states of different areas: 

4/4k) = VJ>1, (70) 

V e v e 

A=pA=\j) = i^+i). 

V e v e 

Nevertheless, we will not use a lot these renormalized operators and we will focus below on dynamical operators built 
from the original operators e, /, / '. 



B. "Black Hole" Creation Operator 



Before studying the dynamics of the spin network states on our 2-vertex graph, we would like to explain further 
the physical interpretation of our homogeneous states \J) oc (P) J |0). In fact, these states are maximally entangled 
states between a and /?. More precisely, let us consider the density matrix corresponding to this pure state and define 
the reduced density matrix on the system a obtained by tracing out the system /3 : 

p^\J)(J\K(P) J \Q)(0\f J , Trp=l, (71) 
Pa = Trap = (j) {J 3)p. 

The main point is that the reduced density matrix on a is the totally mixed state on the Hilbcrt space of iV-valent 
intertwincrs with fixed total area J, 

p a OC I a = I (J) . 

N 

To prove this, we first check that the density matrix p a commutes with the U(7V)-action on Wd' . This comes 
directly from the fact that the full density matrix is by definition invariant under the U(7V)-action acting on both 

vertices, p = (U <S>U) p (U^ 1 (8 U^ 1 ) as matrices on <£> Ti^ . Taking this into account, we have : 

MU G U(iV), UpaU- 1 = Trp (U ® I) ( o(f/" 1 ® I) = Tr^ (I ® U^pil ® U) = Trp p = p a , (72) 

since the trace on "H^f ^ is obviously invariant under unitary since it carries an (irreducible) representation of U(iV). 
Once we have checked that p a commutes with unitary transformations, we can conclude by Schiir lemma that it 
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is necessarily proportional to the identity on the intertwiner space since this space carries (once again) an 

irreducible representation of U(iV). Thus, taking care of the normalization, we conclude that: 

Pa = 7^— V J) ' ( 73 ) 

J->N,J 

This means that the we have a mixed state with maximal entropy on the system a and that our initial pure state is 
maximally entangled. 

In previous work [loj . it was shown that the entropy of this state, given by the logarithm of the dimension of 

the intertwiner space fiffi, grows linearly with the total boundary area J in the large scale limit, and thus behaves 
holographically. It is tempting to call this a "black hole" state in our very simplistic framework, since we have a pure 
state on the coupled system with a "maximal ignorance" state for each of the two sub-systems. From this perspective, 
we can dub our operators (P) J as "black hole" creation operators. 

To insist on this point of view, we remind that it has been recently argued from various perspectives that the Hilbcrt 
space of states for a (quantum) black hole in loop quantum gravity is the space of intertwiners with the same fixed 
total area [l4T - fl8| . This supports our claim here that the induced reduced density matrix p a by the pure homogeneous 
state | J) represents a physical configuration similar to a black hole. 



C. An Ansatz for Dynamics 

Now, we would like to study the dynamics on this 2-vertex graph and propose the simplest U(./V)-invariant ansatz 
for a Hamiltonian operator: 

H = Xe + (af + ap). (74) 

As explained above, the operator e does not affect the total boundary area, [E, e] = 0, while the operators / and /t 
respectively shrink and increase this area, [E,f] = —2/ and [E, /t] = +2/t. The coupling A is real while a can be 
complex a priori, so that the operator H is Hcrmitian. We can relate this Hamiltonian operation to the action of 
holonomy operators acting on all the loops of the 2-vertex graph using the formula (|42[) . From this point of view, our 
proposal is very similar to the standard ansatz for the quantum Hamiltonian constraint in loop quantum gravity [l[ 
and loop quantum cosmology (e.g. (l9l. [20|V 

This Hamiltonian is quadratic in J since it is of order 4 in the oscillators, and we can give its explicit action on the 
basis states of 2 T-Li nv : 

H\J) = a2 v /j{J + 1){N + J- l)(N + J - 2) | J - 1) 
+A2(J+l)(iV + J-1)\J) 

+<t2^(J + 1)(J + 2)(JV + J)(N + J- 1) |J + 1) . (75) 

We will study the properties, especially its spectral properties, in detail in the next section. This operator is particu- 
larly interesting because it allows a direct link with the framework of loop quantum cosmology, as we will discuss in 
section |PVl In the following, we will restrict ourselves to the case of real couplings, a £ R. 

Looking for the eigenvalues and eigenvectors of this Hamiltonian, H J2j a j\J) = 12j a j\J}' ^ s translated into 
a second order recursion relation 5 for the coefficients aj : 

VJ>1, ay/ {J + 1)( J + 2){N + J){N + J - 1) a J+1 +A [(J + 1)(JV + J - 1) - 0\ aj 

+ay/J(J + 1)(N + J — 1)(N + J - 2) = 0. (76) 

We will study in detail the solutions to this equation and the resulting spectral properties of the operator H in the 
next section fill Dl Nevertheless, we can already make some simple comments about it. First, although it is a second 



5 This second order recursion relation is very similar to the recursion relation for the {6j}-symbol of the recoupling theory for SU(2) 
[2lL I *2 2f _ Such a relation can almost be expected since both our states \J) and the {6j'}-symbol define distributions on the space of 
SU(2) intertwiners. Furthermore, directly comparing by hand the two recursion relations, it is possible to have an almost perfect match 
between the coefficients of the recursion relations by setting N = 1 and considering isosceles tetrahedra for the {6j'}-symbol. From this 
relation, we expect three regimes for the coefficients aj corresponding to the three asymptotic regimes of the {6j}-symbol: an oscillatory 
phase (corresponding to true Euclidean tetrahedra), an exponential phase (which can be mapped to tetrahedra embedded in the 2+ld 
Minkowskian space) and a critical regime (which corresponds mathematically to tetrahedra with vanishing squared- volume) . 
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order recursion relation, the initial equation at J = determines the coefficient u\ in terms of u§. Thus we only need 
one initial condition cto and the space of solutions to this recursion relation is one-dimensional once /3 is fixed. The 
second remark is that since it is a second order recursion relation, we expect the coefficients ctj to be oscillatory, at 
least for large J's. Considering the leading order in J in the asymtptotics, it is fairly easy to see that the coefficients 
of the recursion relation only depends on the couplings A and a so the frequency of the oscillations does not depend 
on N and the eigenvalue /? but is expressed only in terms of A and a : 

cosuj = — A/2<7. (77) 

This means that we have a priori three dynamical regimes, depending on the value of the couplings. Taking A > for 
the sake of simplicity, we distinguish the three following cases: 

1. the oscillatory regime when \a\ > A/2 : the frequency uj is real and we expect the coefficients aj to oscillate as 
cxp(iwJ) as J grows large. 

2. the discrete regime when |er| < A/2 : the frequency u is purely imaginary and the coefficients a,j to diverge as 
J grows to infinity. Actually, they always diverge except for a discrete set of choices of eigenvalues (3. Then the 
spectrum of H turns out to be discrete and positive. 

3. the critical regime when a = ±A/2 : the frequency vanishes when a = —A/2 and is equal to it when a = +A/2. 
In the case a = —A/2, our model is very similar to the flat FRW model of loop quantum cosmology (e.g. [23| |) 
and we expect similar asymptotics for our eigenstates. In particular, we will show that the spectrum of H is 
continuous, real, and positive. 

We will study the details of these three regimes in the next section. Moreover, we will discuss later, in section HVBl 
the relation between our couplings A and a and the cosmological constant, comparing these three regimes respectively 
to the three cases A < 0, A > and finally A = 0. 

The ansatz given above is the most general U(N) invariant Hamiltonian (allowing only elementary changes in the 
total area), up to a rcnormalization by a independent factor. Therefore, we can also propose rcnormalizcd Hamiltonian 
operators based on the renormalized operators considered in the previous section. For instance, we define both: 

h= - 1 H — 1 = = XZ + (crA_ + aX + ) e sl 2 , (78) 

y/E + 2{N-l) y/E + 2(N-l) 



h= -=H—==\l + -={af + ap)-=. (79) 
veve ve ve 

Let us underline the fact that we study the action of these Hamiltonian operators H, h and h on the U(A)-invariant 
space 2 'Kinv l nevertheless they are generally well-defined on the whole space of spin network states 2 H.. Let us also 
remind the fact that the parameter N giving the number of edges of the graph has completely disappeared from 
the action of these renormalized Hamiltonians h and h: the dynamics on the homogeneous sector does not depend 
mathematically on N. 

In the context of loop quantum gravity or loop quantum cosmology, these various (re)normalization for the Hamil- 
tonian operator would correspond to different densitized versions of the Hamiltonian operator. Comparing to loop 
quantum cosmology e.g. [liij . our initial Hamiltonian operator seems to correspond to the evolution operator while 
our renormalized Hamiltonian h corresponds to the gravitational part of the LQC Hamiltonian constraint C glav . 

The main characteristic of h is that its coefficients are linear in the variable J at leading order and, more importantly, 
it is an element in the Lie algebra sZ 2 . Thus we know its spectral properties from the representation theory of sZ 2 . More 
precisely, for a real couplings A, a £ R, the generator h = XZ + (aX- + aX + ) lies in the Lie algebra si (2, R) ~ su(l, 1). 
The (Lorentzian) norm of this 3-vector, 

||h|| 2 = A 2 -4| ( t| 2 , (80) 

determines its type. There are three cases : 

1. If A 2 > 4 1 cr | 2 , h is time-like, i.e it generates a space rotation in SU(1, 1) and it is conjugated to Z by a SU(1, 1) 
transformation (up to a numerical factor). In this case, the spectrum of h is discrete. More precisely, since the 
spectrum of Z is N*, then the spectrum of h is sign(A) \ \h\ \ N*. 

2. If A 2 < 4|ct| 2 , h is space-like, i.e it generates a boost in SU(1, 1), its spectrum is continuous and covers the whole 
real line. 
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3. If A 2 = 4 1 cr | 2 , h is null-like and its spectrum is also continuous and of the same sign than A. 

For more details on the spectral analysis of su(l, 1) algebra elements and their eigenvectors, the interested reader 
can refer to (24| . The important point to underline here is that the three coupling regimes for the renormalized 
Hamiltonian h are exactly the same as for the original Hamiltonian H . This is very similar to the interplay between 
the two operators in LQC, the evolution operator and the gravitational contribution to the Hamiltonian constraint 
Cgrav (see e.g. 

The advantage of this Hamiltonian operator h over the initial Hamiltonian H is that its properties are much 
better known. For instance, its spectral properties are well-established and we know its exact discrete spectrum in 
the weak coupling case A 2 > 4|er| 2 . Also the exponentiation of Lie algebra elements is under control and we know 
that h will generate finite SU(1, 1) transformations. Nevertheless, we will also analyze the properties of H since this 
choice of Hamiltonian seems to correspond to the operator used in LQC. 

The main feature of the last operator h is that its coefficients are of order in the area J. Actually, its action is 
truly simple: 

h\J) = \\J) +a\J-l)+a\J + l), h\0) = |0) (81) 

This operator has properties very similar to the holonomy operator x whose action was given in (|45[) . Restricting 
ourselves to the case of a real coupling a £ R for the sake of simplicity, the Hamiltonian h has a bounded spectrum 
and its eigenvectors (in the generalized sense) are given by plane waves (for more details, see the appendix) : 

h\uS) = (A + 2crcosw)|o;), with = V] sin [(J + l)w] | J). (82) 

JGN 

The eigenvalue (3 = A + 2<rcosu; corresponding to the frequency cj is obviously bounded by the extremal values A±2<7. 
This operator is less interesting than H and h, and is not relevant to the relationship between our framework to loop 
quantum cosmology. 



D. Solving the 2- Vertex Model 



In this section, we study in more details the spectral properties of the Hamiltonian H depending on the couplings 
A and a. In particular, we will establish the existence of the three regimes discussed above. To this purpose, we 
will use the techniques used in LQC to study the properties of the gravitational part of the Hamiltonian constraint 
[13, HH, HH- In the following, we will assume that A > for the sake of simplicity. 

We start by checking the positivity of our operator (which we rescale by a factor 2): 



H\J) = \<p{J) \J)+aip{J) \J-l)+aip(J+l) I J+l) , with 



tp(J) = ( J + 1)(N+ J- 1) 

ip{J) = y/J(J + 1)(N + J — 1)(N + J - 2) 



(83) 



We consider the diagonal part of this operator as our free Hamiltonian Ho | J) = (p(J) \J). Actually it is simply our 
operator e up to a factor 2. Then the off-diagonal components triggering the transitions to the | J ± 1) states are 
considered as the interaction: 



H = XH Q + I with I\J) = aip(J)\J -l)+aip{J+l)\J+l). 
Result 1. For a weak coupling, |er| < A/2, the Hamiltonian H is positive, \fv (v\H\v) > . 



(84) 



Proof. Writing the vector v in the fundamental basis, v = ^2jVj |J), we compute explicitly the action of the Hamil- 
tonian: 



MH\v) = 



j 



\ip(J)\vj\ 2 + aip(J)vj-ivj + aip(J + l)vj+wj . 



Then to show the positivity, since tp( J) > and thus the free term is always positive, we just need to show that the 
interaction term is smaller in than this free term. First we use the obvious inequality, < (|wj| 2 + \vj + i\ 2 )/2, 

which means that: 



^ o-il>(J)vj-iVj + a^(J + l)v j+ ivj 



< 



W\ £>(J) + V0/ + i))k/| 



18 



Now, we only have to check that < (ip(J) + tp(J + 1)) < 2<p(J) holds for all J's. This is straightforward to show 
by squaring this inequality. Therefore, we can conclude that, as soon as 2|er| < A, the Hamiltonian is positive. On 
the other hand, it is easy to see that as soon as 2|er| > A, then the expectation value (v\H\v) can take arbitrary real 
values. 

□ 

Result 2. We can bound the interaction term by the free Hamiltonian: 

W, ||/v|| < 2|cr| \\H Q v\\ + \a\ \\v\\ . (85) 
Proof. We follow the same steps as in the previous proof: 



\Iv\\ 2 = 



o-<P{J)vj I J - 1) + ™l>V + \ J + 1) 
J J 

< 2|a| 2 ^(|^(J)| 2 + |^J + l)| 2 )| V .7| 2 , (86) 
J 

where we used that for an arbitrary couple of vectors v, w, we have the norm inequality | |t> + w;| | 2 < 2(||u|| 2 + ||u;|| 2 ). 
Then we compute explicitly: 

4>{J) 2 + ip(J + l) 2 = 2(J + 1)(N + J - 1)(J 2 + NJ + N) = 2(p(J)((p(J) + 1) . 

Thus, inserting this formula back in the previous inequality, we get: 

HHI 2 < 4|a| 2 (||ii' 7;|| 2 + (7;|H 0t ;)) 

< 4H 2 (||ff,HI 2 + IMII|£MI) 

< 4H a f||F «|| + J|H| > ) , (87) 
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which allows us to obtain the bound. 



□ 



Since the operator Hq is obviously essentially self-adjoint, and since H = XHq + /, we can apply the previous 
inequality in the small coupling case, 

M < \ => \\Iv\\ < \\XH v\\ + W\\\v\\, 

and conclude that: 

Result 3. The Hamiltonian H is essentially self-adjoint as soon as 2\o~\ < A. 

This allows to prove the existence of the discrete regime when 2|<r| < A : 

Result 4. When 2\a\ < A, the spectrum of H is discrete and strictly positive. 

Proof. Assuming that 2|er| < A, we consider the following decomposition of the Hamiltonian: 

H=(\-2\a\)H + (2\a\H + l). 

Using the result 1, we know that the operator (2|cr|_ffo + I) is positive. Thus we have the operator inequality 
< (A-2|cr|)ff < H. 

Finally, we use the facts that H is essentially self-adjoint for our choice of coupling and that (A — 2|cr|)J?o has a 
discrete positive spectrum and we apply the lemma in [25| (p-16) to conclude that H has a discrete spectrum. 

We do not repeat the rigorous proof of this lemma here. Nevertheless, the intuitive idea is that the subspace 
generated by eigenvectors of H with eigenvalue inferior to an arbitrary value fJ m is automatically included in the 
subspace of of eigenvectors of (A — 2|<r|)i?o with eigenvalue inferior to the same max value j3 m . Since the latter has a 
finite dimension, the first also does. 

□ 
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Now we turn to the critical regime when 2\a\ = A. For the sake of simplicity, we fix A = 2, then a is on the unit 
circle, a = exp(—i9) with an arbitrary angle 9 £ [0, 2ir]. We will write H c for the Hamiltonian operator in this critical 
regime. We already know from the previous analysis that H c is essentially self-adjoint and strictly positive. We will 
not provide here rigorous proofs about its spectral properties. We will present instead approximative asymptotics 
for its eigenvectors. For a more rigorous approach, one could use techniques similar than the ones used in loop 
quantum cosmology, performing for instance a Fourier transform on a simpler operator which would differ from H c 
by a trace-class operator [2(| . 

To identify the eigenstates, thus solving the equation aj H c \ J) = ctj | J), we have to solve the recursion 

relation: 



VJ e N, 2cp( J)aj + e +i6 ip(J)aj-i + e~ l V(J + l)a J+1 = (3a.j. 



Since this is a second order recursion relation, we expect a two-fold degeneracy for the eigenvectors. Moreover, following 
the results obtained in the framework of loop quantum cosmology, we expect an asymtptotics of the eigenvectors of 
the type aj ~ exp(ifc In J)/ \[~J [23|, H(| [27j • This is compatible with the intuition that the frequency to of oscillations 
in J vanishes in the critical regime, since we now have oscillations in In J which are much slower. Nevertheless, due 
to the initial relation for J = 0, the coefficient a\ is entirely determined by the initial coefficient ag. Therefore, the 
whole sequence (ctj)j is determined by a single initial condition a.Q. This kills the two-fold degeneracy and implies 
that the eigenvalues have a trivial multiplicity of 1. 

Following this line of thought, let us insert the following ansatz in the recursion relation: 



otj 



(-1)- 



g2# J ^ik In J 



(89) 



with k € R. We insert this ansatz in the relation above (|88j) and we study the large J behavior pushing the 
approximation to second order in J, that is in 0(1) : 



cp(J) 



otj-i 



CiJ+1 
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(N - l)J 
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A straightforward calculation allows to check that this ansatz does indeed satisfy approximatively (at the order 0(1)) 
the recursion relation for an eigenvalue easily expressed in terms of the frequency k : 



P = t + k*. 



(90) 



As expected, this leads to a strictly positive spectrum, with a two- fold degeneracy since (3 is in k 2 . This gives us a 
first clue about the spectral properties of H c . 

Conjecture 1. Considering the critical regime with 2\a\ = A with A > 0, then the (generalized) spectrum of our 
Hamiltonian H is continuous and covers ]-|,+oo[. All the eigenvalues have a multiplicity 1. 

We can check the orthogonality of the eigenvector ansatz : 

J(k-k)\nJ ^ I dx X e i( - k ~ k ^ lnX ' 



dye i{k - k)v (x 5{k-%). 



(91) 



Now, in order to get a precise prediction for the asymptotics of the eigenvectors, we need to check how the initial 
condition collapses the apparent two-fold degeneracy to a trivial multiplicity. A priori, it does not make sense to 
apply the initial condition relation between ao and <x\ to our asymptotic ansatz for large J. Moreover, our ansatz 
clearly diverges for J = 0. 
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We set the initial condition ao to be real, for instance we choose ao = 1. A first remark is that all the coefficients 
a,j are then real too up to the exp(idJ) phase due to the coupling a since the coefficients of the recursion relation 
and the eigenvalue /3 are themselves real. Therefore, a more precise ansatz for the asymptotics of the eigenvectors of 
large areas J is: 

(-lV 

aj ~ A(k) ±-;L- e i9J cos(fc In J + cj)(k)), (92) 
v J 

where (j>(k) is a phase depending on the eigenvalue and to be determined, and A{k) is the amplitude of the oscillations 
which depends only on k and is proportional to the initial condition a - Unfortunately, we do not yet have any 
analytical predictions or simple numerical fits for A{k) and 4>(k). Nevertheless, we provide a couple of plots to 
illustrate the asymptotical behavior of the eigenvectors in figiH 




FIG. 3: Plots in the critical regime for the coefficients aj for values N = 4, er = — 1 and two different frequencies k = 3 on 
the top and k — 15 below. The plots of the left are computed exactly using the recursion relation, while the plots of the right 
are given by the asymptotic formula. We notice that we have the right oscillatory behavior and the correct overall scaling in 
l/y/7, but we still miss approximative predictions for the full amplitude and the initial off-shift of the oscillations. 



We now turn to the strong coupling case, with |A| < 2\a\ which correspond to our oscillatory regime. The recursion 
relation corresponding to the eigenvalue equation, H J^j a J \J) — ft a J \J) is as before: 

VJ e N, \<p( J)aj + aij){ J)aj_i + a^( J + l)a J+1 = potj. (93) 

As pointed out earlier, in this strong coupling, with |A| < 2\a\, we can define a real frequency u in terms of the 
couplings: 

A + ae-^ + ae +tu = 0. (94) 
This allows us to pick up an oscillatory ansatz for a leading order approximation of the eigenvectors: 

aj = -^e- (95) 

Inserting this in the recursion relation and pushing the approximation to the next-to-next-to leading order in J, that is 
to the order 0(1), we use the previous approximations for the factors tp( J) and ip{J), and the following approximations 
for the coefficients: 



aj±i ~ aje ±m I 1 =F ~ 



1 1 ± c 



J J 2 

A straightforward calculation gives us the following eigenvalue: 

P = (f -c) (ae-^-ae+n- (96) 

This is a surprising result: /3 is complex for complex values of the parameter c, and even purely imaginary when ceK. 
Moreover the coefficients aj are I? since they go as 1/J for large J's. Therefore, we have normalizable cigenstatcs 
corresponding to non-real eigenvalues. This leads us to the following statement. 
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Conjecture 2. Considering the oscillatory regime with 2|er| > |A| then the Hamiltonian H is not essentially self- 
adjoint and its spectrum covers the whole complex plane C. 

This claim is very similar to what happens to the Hamiltonian constraint in loop quantum cosmology for a positive 
cosmological constant A > [l9| . 

To conclude, the spectral properties of the Hamiltonian H are well-understood, although we do not have rigorous 
proofs in the critical and oscillatory regimes and we do not have the exact discrete spectrum in the weak coupling 
regime. Nevertheless, we know the spectrum of H in all three cases and we have the explicit leading order asymptotics 
of the eigenvectors in the critical and oscillatory regimes. 

However, at the end of the day, it might still be better to work with the s^-Hamiltonian h whose spectral properties 
are entirely known and proved. For instance, its discrete spectrum in the weak coupling is explicitly known. Moreover, 
we know how to exponentiate it into SU(1, 1) group elements, which would define the finite time evolution of the 
system. 

IV. COMPARING WITH LOOP QUANTUM COSMOLOGY 

As already shown in @, 0] , the 2-vertex graph is a perfect setting to derive a quantum cosmology sector from the full 
loop quantum gravity theory. More precisely, the authors of these papers show how their model was related to loop 
quantum cosmology. In this section, we similarly underline the explicit relation with the loop quantum cosmology 
framework. 



A. Almost the Same Mathematics 

Comparing our Hamiltonian operator(s) to the loop quantum cosmology framework (e.g. we immediately 

notice some substantial differences: 

• The basis states are labeled by a volume variable v instead of our total area variable J, and the (modified) 
holonomy operator induce constant shifts in the volume v — > v ± 4 which can not be interpreted as constant 
shifts in the area assuming that the area goes as usual as J ~ v 3. 

• The volume variable is a priori real, v G R, although supcrselection rules restrict it to a discrete lattice. 

• The volume v is a priori not bounded by below and is not restricted kinematically to be positive. The restriction 
to a v > sector occurs at the dynamical level. 

From this perspective, our framework is much closer to the original loop quantum cosmology setting initially intro- 
duced by Bojowald [28j, |29j]. Despite these differences, the mathematics of the gravitational part of the Hamiltonian 
constraint in loop quantum cosmology (LQC) are very similar to the mathematics of our Hamiltonian operator. 

Indeed, there exist various orderings which have been investigated in isotropic LQC, for example the APS prescrip- 
tion [3(| or the sLQC simplified model ("s" standing for solvable) [3l| or the MMO ordering (32[, but they all have the 
same leading order structure at large volume. More precisely, looking at the flat case with a vanishing cosmological 
constant A = 0, the LQC framework usually studies two related operators, the gravitational part of the Hamiltonian 
constraint C gr and the evolution operator 0, which differ from one another by a density factor. Considering the 
asymptotic regime, v 1, the action of these two operators reads as (e.g. (27j): 

C gr \v) ~ (A(v + 2) + A(v - 2)) \v) - A(v + 2) \v + 4) - A(v - 2) \v - 4), (97) 

= VvC gr y/v => 9 \v) ~ v(A(v+2)+A(v-2)) \v)-y/v(v + 4) A(v+2) \v+A)-^v(v - 4) A{v-2) \v-A) , (98) 

with A(v) oc v at leading order in the asymptotic limit. Comparing these expressions to the action of our Hamiltonian 
operator(s) |83|) . it is obvious that they have the same action at leading order if we substitute the LQC variable v by 
our total area label J, or more exactly by setting v = 4 J. Indeed the coefficients of C gr grow as v and the coefficients 
of 9 grow as v 2 : 



C gr \v) (x 2v\v) — v\v + A)-v\v — A), 
Q\v) cx 2v 2 \v) - v 2 \v + 4) - v 2 \v - 4) , 



(99) 
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up to an overall factor. This leads us to two conclusions. First, the evolution operator O corresponds to our 
Hamiltonian operator H, whose coefficients grow as J 2 , while the LQC Hamiltonian constraint corresponds to our 
s/2-Hamiltonian h, whose coefficients grow as J at leading order. Therefore, the spectral properties of the LQC 
operators and of our Hamiltonian operators will be very similar and we can apply the techniques already developed 
in LQC to our framework (2(J [25|, [26[ , and vice- versa. 

Second, these LQC operators for the flat case A = clearly correspond to our critical regime with a = —A/2. One 
can directly compare the properties of the LQC operators to ours in this regime and we see, for example, that the 
eigenstates have exactly the same leading order asymptotics respectively in exp(ifc In v)y/v and exp(ifc In J) \/~J . 

Finally, this identification of the vanishing cosmological constant case A = with our critical regime leads to 
the natural idea that the various regimes of our operators correspond to the different signs of A. We discuss this 
correspondence in more details below. 

At the end of the day, the present comparison between LQC and our framework is rather superficial. A deeper un- 
derstanding of the precise relation at the physical level could be achieved when matter or anisotropies/inhomogeneities 
are included into our 2-vertex model. In particular, this would allow to determine whether our choice of area vari- 
able versus the standard LQC choice of a volume variable still leads to sensible physical predictions. This will be 
investigated in future work [351 ] . 



B. The Physical Meaning of the Couplings and the Cosmological Constant 

We would like to look at the way the cosmological constant is usually included in isotropic loop quantum cosmology. 
Following [l9| , the cosmological constant simply leads to a diagonal volume contribution to the Hamiltonian constraint: 

C gr \v) = {A{v + 2) + A(v - 2)) \v) - A(v + 2)\v + 4} - A(v - 2)\v - 4) - AV \v) , (100) 

where we have taken v > for the sake of simplicity. The volume operator V is obviously assumed to be diagonal: 

V\v) = V v\v), 

and the coefficient A(v) is assumed to be linear in v in the asymptotic limit: 

A(v) ~ 2A v. 

The units Aq and Vb ar( 3 usually defined in terms of the gravitational constant G, the Planck constant ft, and the 
Immirzi parameter 7. The interested reader can find the explicit expressions in [l9j . Since the basic shifts in the label 
V are ±4, it is natural to identify v = 4J in order to write explicitly the correspondence with our U(JV) framework. 
Moreover, we study 8 = ^JvC gr ^Jv instead of C gr . Translating it in terms of the variable J, we get at leading order : 

S\J) ~ 16(4A - A Vb) J 2 I J) - 32A a (J + ^)y/J(J + 1) \J + 1) - 32A (J - ^)y/J(J — 1) | J — 1). (101) 

We compare it directly with the expression of our Hamiltonian operator H given in (|83j). The two expressions match 
exactly at leading order in J. It is now clear that when A = 0, we are in the critical regime. On the other hand, as 
soon as A ^ 0, we leave the critical regime and enter either the discrete or oscillatory regimes. More precisely, we 
have the following correspondence between our couplings and the LQC parameters: 

A = 16(4A - AV ), a = a = -32A . (102) 

Thus when A = 0, we have a = —A/2 and we are in the critical regime. When A > 0, but still close to 0, we have 
< A < 2\a\ and we are in the strong coupling regime with a non-essentially self-adjoint Hamiltonian operator H. 
Finally, when A < 0, we have A > 2|er| and we are in the weak coupling regime with a strictly positive discrete 
spectrum of the Hamiltonian H. These 3 cases perfectly fit with the results obtained in LQC fl9l |33| . 
This correspondence is however very rough. It needs to be refined on at least two points: 

• Here we made the substitution v — > 4J at the purely mathematical level. Of course, physically, we can not 
replace the volume variable v by the area label J. Thus we would need to understand better the role and the 
validity of such a substitution and in particularly we should study the expectation value of the volume of our 
isotropic states | J) . 

• We need to include matter in our 2-vertex model to truly compare it with the full Hamiltonian constraint and 
the exact LQC results. 
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Besides these set-backs, it is clear that the cosmological constant contributes a diagonal term to the Hamiltonian in 
the | J) basis (since the value of the volume should be given uniquely by the total area J for isotropic configurations) 
so that shifting A is naturally reflected in the relative values of the couplings A and a. From there, it is natural 
to identify the three regimes of our Hamiltonian, discrete or critical or oscillatory, respectively to the three sectors 
A <=> 0. 

From this perspective, we have a physical interpretation to arbitrary real values of the coupling a. We would still 
need to understand the physical meaning of the possible complex phase of that coupling. 

C. Going further with LQC? 

So far, we have shown how our framework resembles the LQC framework and how the LQC techniques can be 
used to study the properties of our Hamiltonian. Moreover, the LQC framework is much more developed than ours 
and includes matter fields and investigates various cosmological settings. Thus, a natural question is: what can our 
approach bring to LQC? Wc would like to present two points on which our approach can bring something new to 
LQC. 

Our first point is technical. It's the remark that our operators acting on the homogeneous subspacc form a 5I2 
algebra. This allows us to characterize our homogeneous sector as an irreducible unitary representation of SL(2,K) 
and to construct a Hamiltonian operator h as a s^-Lie algebra element whose spectrum and eigenvectors are explicitly 
known (in all cases of the couplings). Since the structure of the actions of our Hamiltonian operator (s) and of the 
LQC Hamiltonian constraint are very similar, wc can expect that a particular choice of ordering in LQC would also 
lead to a sZ2-structure. This would allow for a more refined algebraic characterization of the LQC Hilbert space and 
of the spectrum of the Hamiltonian of isotropic LQC. 

Our second point has deeper consequences. It's on the crucial role of U(A r ) as the symmetry at the quantum level 
reducing our full Hilbert space of spin network states to the subspace of isotropic and homogeneous states. Indeed, the 
standard loop quantum cosmology framework relies on the implementation of the requirement of isotropy/homogeneity 
at the classical level: we reduce the degrees of freedom at the classical level, and then quantize the reduced phase 
space using loop quantum gravity techniques. Here the procedure is done the other way round: we start with the full 
loop quantum gravity and reduce the degrees of freedom directly at the quantum level by first selecting the underlying 
2-vertex graph with an arbitrary number N of edges and then imposing the invariance under U(AT) in order to reduce 
our Hilbert space to isotropic and homogeneous states. From this viewpoint, the U(AT) symmetry could be the tool, 
missing up to now, necessary to perform the symmetry reduction directly at the quantum level from loop quantum 
gravity to isotropic cosmological situations. 

V. GOING BEYOND THE 2-VERTEX MODEL 

In this section we will comment on the main limitations of our model and some possible generalizations of it. More 
precisely, we will propose a tentative model with 3+A^ vertices where, in principle, it could be possible to study certain 
additional physically relevant processes. 

Let us start by sketching, in first place, some possible generalizations of our model (in order of increasing difficulties): 

• There are neither inhomogencities nor anisotropics in our model since we have restricted ourselves to work on 
the U(A r )-invariant space, which restricts to isotropic states (| J)) only labeled by the total boundary area J. We 
could relax this requirement and study the dynamics of non-U(A^)-invariant states. Since our dynamics is \J(N)- 
invariant, irreducible representations of the U(N) group would still be stable under the time evolution. It would 
be interesting to check what type of anisotropics and inhomogeneities correspond to each of these irreducible 
U(AT) representations in the spin network space of our 2-vertex graph, use U(iV") to classify anisotropics and 
compare to the standard expansion in terms of angular momentum. 

• We could go further and break our requirement of U(AT) invariance of the dynamics. We could do this by hand, 
by choosing operators which are not invariant, but it does not seem physically motivated to break a symmetry 
for free with no clear purpose. Moreover, the dynamics of the two vertices is exactly the same due to our choice 
of dynamical operators, which are symmetric, so there is no freedom for local fluctuations. Furthermore, there 
is no matter in our simple model. We could add a scalar field on the 2-vertex graph (like in 0). Matter could 
provide us with a source to break the U(iV) invariance and introduce local fluctuations in the model. It would 
also allow to obtain a true (though still very simple) cosmological model and compare it with the results in loop 
quantum cosmology. 
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• Our two vertices are strongly glued with the matching conditions. They grow, shrink and rotate, all together with 
no freedom for different evolutions for the two vertices. We cannot relax this without changing the underlying 
graph and using a more complicated spin network setting. This would allow us to represent physical processes 
such as the rotation of one vertex with respect to another or the evaporation of a (horizon) boundary. Here we 
propose to consider the next simplest graph, a 3+iV vertex graph, which allows to decouple the two original 
vertices and recover the independent U(iV)-actions on the intertwiner spaces of both vertices. 

We postpone these points to future investigation. Nevertheless, we would like to describe here the 3+N vertex 
graph that we just mentioned. It is the simplest graph beyond the 2-vertex graph and we will sketch how it would 
allow a richer physical arena in which more processes could take place. 

The main physical limitations of the 2-vertex graph considered so far are direct consequences of the gluing between 
both vertices. We have seen how this model is good enough to describe homogeneous cosmological situations, and 
it would be also possible to consider the introduction of inhomogeneities by relaxing the U(iV) symmetry, as just 
commented above. However, an interesting point of view is to regard the two vertices as representing an "in" and an 
"out" regions. Adopting this point of view, we would like to be able to describe physical situations in which the "in" 
region has some kind of dynamical evolution with respect to the "out" region. Examples of these processes could be 
evaporation of the "in" region towards the exterior, or rotations of the interior. The present 2-vertex model turns out 
to be too simple to describe this kind of physical situations. For instance, in particular, the "out" region cannot exist 
without the "in" region, so in an evaporation process (e.g. as in the case of a physical black hole) causing the "in" 
region to shrink and eventually disappear, the "out" region would also be forced to disappear. These limitations are 
imposed by the graph itself. Thus in order to introduce enough freedom to study this kind of processes, we propose 
to consider a slightly more sophisticated graph. 

This generalization consists in adding a third vertex to the graph, which will be part of the "out" region. This new 
vertex 17 is linked to the previous graph by N edges, each one meeting one of the old edges at a new trivalent vertex 
<Zj. The resulting picture is shown in figure [4j The interior region continues to be the vertex a and the boundary 
is the surface surrounding it. The N trivalent vertices are considered as part of the exterior region, but can also be 
interpreted as representing the physical boundary between the "in" and "out" regions. 




FIG. 4: The 3-vertex graph: the three vertices a, ft, and Q, are linked by N edges and we have TV additional (auxiliary) trivalent 
vertices. 

We will consider operators acting on the new vertices in analogy to the ones we already have for the simpler 2-vertex 
model. Thus, we have operators E^' and acting on the vertex fi as well as Eh. and i* 1 }? , with i = 1, . . . , N 
and j, k = a, p, f2, acting on the N trivalent vertices. 

Now, thanks to this new vertex fi it is possible to shrink the area surrounding the "in" region without resulting in 
a shrinking of the "out" region too. One can even think of the "in" region totally disappearing, and the "out" region 
would still make sense. Thus, in this framework it is possible to introduce radiation-like operators. As a first attempt 
we propose an operator of the kind 

R^=Fl?E^F^E^. (103) 

The action of this operator over the graph is to decrease the spin associated to the edges i and j of the interior 
vertex, and to increase the spins of the same edges of the vertex Q. In this sense, the action of the operator decreases 
the area of the boundary by "radiating" some spin units to the exterior region. The vertex f2 could be interpreted in 
this framework as carrying the radiated field. 

In a similar way, one can introduce operators generating rotation-like transformations, which would change the 
representations living on the N legs of the vertex a but without changing the total area of the boundary around a 
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FIG. 5: The action of the radiation operator on the left and the rotation-like operator on the right: we shift the spins on the 
legs i, j around the vertices a and fl without affecting the spins living on the legs around the vertex /3. 

and without affecting the representations living on the N legs of the vertex /3. It will be necessary to check whether 
the commutators of these radiation operators and rotation-like operators satisfy a closed algebra. We postpone the 
study of the details of these generalizations, as well as the implementation of the dynamics for this new 3+iV model, 
to future work [34j . 

VI. CONCLUSIONS AND OUTLOOK 

The U(7V) framework recently introduced in [ij-flil] represents a new and refreshing way to tackle the issue of 
dynamics within the loop quantum gravity framework. Our motivation is not to solve the full set of constraints but 
to propose a Hamiltonian operator which describes the dynamical behavior of a given system, following the usual 
procedure even in general relativity. Here, we have presented a simple 2-vertex model -based literally on graphs with 
2 vertices and an arbitrary number ./V of edges- and we have proposed a tentative Hamiltonian for the dynamics of 
this specific system. 

In those previous papers [ij-[Tl|, it was shown that the space of SU(2) intertwiners for loop quantum gravity carries 
a natural action of the unitary group U(A) : the Hilbert space of intertwiners for a given number N of edges and with 
a fixed total area J is identified as an irreducible U(7V) representation. This construction is based on the Schwinger 
representation of SU(2) in terms of harmonic oscillators. Using this tool, a set of operators invariant under global 
SU(2) transformations have been identified : the u(7V) generators and annihilation/creation operators Fij,F^, 
which endow the full space of ./V-valent intertwiners with a Fock space structure. 

In this work we have described a very simple model for loop quantum gravity, based on the 2-vcrtcx graph described 
above. We have the space of intertwiners for each vertex, plus some matching conditions coupling the two intertwiners 
by imposing that each edge i of the graph carries a unique spin label ji as seen from both vertices (these conditions 
reflect that one edge carries a single SU(2) irreducible representation). Taking this into account, we have introduced 
operators e^j , fij , fjj that commute with the matching conditions and constitute the building blocks for the dynamics 
on our 2-vertex graph. Going further, we have identified these matching conditions as the Cartan subalgebra of a 
larger u(N) algebra acting on both vertices, and we have shown that the invariance under this global U(7V) symmetry 
implies the restriction to isotropic homogeneous states | J) which are not sensitive to area-preserving deformations of 
the boundary between a and (3. This provides our U(iV) symmetry with a physical interpretation. This is actually 
one of the key results presented here. Indeed it could be viewed as a key step towards a full understanding between 
the general loop quantum gravity (LQG) framework and the symmetry-reduced loop quantum cosmology (LQC). 
More precisely, the standard framework of loop quantum cosmology is based on a symmetry reduction at the classical 
level followed by a quantization using loop gravity techniques. Since the symmetry reduction is done before the 
quantization, we have lost the connection between the Hilbert spaces of LQG and LQC and we don't know how to 
perform the symmetry reduction directly at the level of the full quantum theory A first hint towards solving this 
issue was provided by the approach of Rovelli-Vidotto which is also based on the 2-vertex graph with N = 4 
edges. They propose to consider the total area J = Y^iJi as the homogeneous degree of freedom and to interpret 
the individual fluctuations of the spins ji on each edge as the anisotropy /inhomogeneous degrees of freedom. They 
implemented this separation of variables by hand when studying the dynamics of their system. Here we went one 
step further and identified explicitly the symmetry which selects the homogeneous states. Now the reduction to the 
Hilbert space of isotropic homogeneous states can be achieved through a simple U(A) group averaging. We hope that 
this technique will be relevant to understanding how to derive quantum cosmology from the LQG framework directly 
at the quantum level. This will be investigated in more details in future work [35j. 

Another interesting point that we would like to stress is that we have studied the relation between our operators 
and the usual holonomy and flux (derivation/grasping) operators standardly used in LQG. This provides an explicit 
bridge between our formulation and the standard way to formulate dynamics in LQG. From this point of view, we 
hope that our formalism will be interesting to study the general issue of the LQG dynamics on arbitrary graphs, more 



2G 



complicated than the 2- vertex graph used in the present work. 

The model that we have developed here has a simple enough algebraic structure to allow us to construct a Hamil- 
tonian that respects isotropy and evolves homogeneous states into homogeneous states. In order to be able to build 
such a Hamiltonian, we have defined U(A)-invariant operators e and / as linear combinations (summing over all the 
edges in order to have homogeneity) of our building block operators Sij and . It turns out that these invariant 
operators simply satisfy a sZ2-algebra. Using the operators / and /' as (area) annihilation and creation operators, we 
have studied their action on an homogeneous state and we have discussed the interpretation of the creation operators 
(f^) J as "black hole" creation operators. This allowed us to define a U(7V)-invariant Hamiltonian, which is unique 
up to a rcnormalization by a factor depending on the total boundary area. Here, we have focused on two precise 
choices of normalization, which allow a direct comparison with the Hamiltonian constraints used in LQC, and we 
have studied their spectral properties. We hope to be able to apply in the future this analysis of their properties to 
further develop the derivation of LQC from our 2-vertex LQG model [35j . following the path opened by Rovelli and 
Vidotto H0j. 

To conclude, although we have been dealing with a simplified model, and one has to be cautious in order to 
extrapolate the results to more general frameworks, in our opinion, we have obtained sensible results and it seems 
that the connections of our 2-vertex quantum model with LQC and LQG arc fully justified. From this perspective, we 
believe that it is plausible to use this kind of simple models as a testbed for certain dynamical problems within the LQG 
and LQC theories. A possible extension of the work would be the implementation of anisotropies/inhomogeneities 
and the study of the stability of homogeneous states under the evolution given by our Hamiltonian. Finally, we also 
discussed the possibility of working on a more complex graph and switching to a 3 + iV-vertex graph, which would 
allow to construct operators representing physical rotation and evaporation and thus to simulate the dynamics of 
physical quantum black holes in LQG. 
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Appendix A: Diagonalizing the renormalized Hamiltonian 

In this appendix, we study the action of the renormalized Hamiltonian operator h (with a £ K), whose action on 
the homogeneous space 2 Tii nv is given by : 

h\J) = A|J) +cr|J- 1) +o-|J+ 1>, h\0) = |0) 

Looking for eigenvectors h J2j a J \J) = ft J2j a J wc obtain a simple second order recursion relation on the 
coefficients a,j of the corresponding state: 

Aa + ercti = (3a , (Al) 
Act j + era j + craj+i = /3a j , VJ > 1 . (A2) 

Putting aside the initial condition equation (|A1[) and focusing on the recursion relation (|A2[) . we see that the plane 
waves a,j = cxp(iuJ) with arbitrary frequency satisfy this second order recursion : 

\e luJ + ae iu( - J -V + m"'^ 11 = (A + 2crcoso;) e luJ . (A3) 

The eigenvalues (3 = (A + 2ctcosw) are obviously bounded and the spectrum is doubly degenerate since the two 
conjugate plane waves exp(i ± ujJ) correspond to the same eigenvalue. Now, we need to take into account the initial 
condition (|A1[) and this leads to a non-degenerate spectrum with a single plane wave for each eigenvalue/frequency. 
Indeed, considering the ansatz aj = a cos a; J + b sin u> J, the initial condition equation reads : 

Xa + <t (a cos lj + b sin uj) = f3a = (A + 2c cos w) a. (A4) 
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This relates the two coefficients a and b of the plane wave, 

fesinw = acosw, (A5) 
which leads to a unique eigenvector for each frequency uj: 

h\oj) = (X + 2a cos lo) \lo) , with \oj) = (cos ui sin w J + sin uj cos w J) | J) (A6) 

= I>in[(j + iH l J ) = sinw ]L^HI J >- 

Jew ./ 

Due to the invariance under parity uj — > — w, we restrict ourselves to positive frequencies uj £ [0,7r]. Then it is 
straightforward to check that these eigenstates are orthogonal to each other: 

(lu\u)) = sin [( J + l)uj] sin [( J + 1)uj] = sin nuj sin nuj = — S(uj — uj). (A7) 
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